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^ ' Abstract 

Two electrons move in a quasi one-dimensional wire under the influ- 
ence of a short-range interaction. We restrict Hilbert space to those states 
^ ' where the two electrons are close to each other. Using supersymmetry, we 

' present a complete analytical solution to this problem. The two-body in- 

teraction affects the density of states and, thereby, the localization length. 
We derive a criterion for the onset of changes of the localization length 
due to the two-body interaction. 
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1 Introduction 



The eigenfunctions of electrons moving in a random disorder potential in one or 
two dimensions display localization. The influence of the Coulomb interaction 
(or of any other two-body interaction) on localization properties has been an 
open problem for a long time. In 1994, Shepelyansky [Q investigated the motion 
of two interacting electrons in a disordered one-dimensional wire numerically. 
He found that the two-body interaction may strongly increase the localization 
length, independently of the sign of the interaction. This surprising result was 
rederived by ImryjS who used Thouless' block scaling picture. Further numer- 
ical work [|[ |[ ^, ^ ^, H H in one dimension aimed at a better understanding 
and quantification of these findings, and at settling a controversy which had 
arisen. As a result, the effect is now well established. In addition. Song and von 
Oppen 1^ have demonstrated that it depends strongly on the distance between 
the two electrons and is strongest when these remain close together. 

It is obviously desirable to extend these investigations to two-dimensional 
systems. At present, numerical approaches do not seem capable of handling 
the ensuing difficulties. However, the work of Refs. |[ |j yields an analytical 
expression for the change of the localization length ({U) under the influence of 
a two-body interaction with typical matrix element U. The result 
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with B the bandwidth due to the disorder potential, is independent of the sign 
of the two-body interaction and is claimed to hold both in one and in two di- 
mensions. Numerical results in one dimension although in line with the 
sign-independence, do not fully support the dependence predicted analyti- 
cally. Moreover, the effect is numerically found to be strongest when the two 
electrons move at short distance a feature which is also not covered by the 
analytical expression. Finally, the result (|l|) has either been derived heuristi- 
cally iQ, or it is based on statistical assumptions [D which, although plausible, 
are not obviously valid. 

This situation calls for a novel approach to the problem. In the present pa- 
per, we extend the study of localization properties of two interacting electrons 
to the case where the electrons move in a quasi one-dimensional disordered wire. 
This case is obviously more realistic (and numerically much more difficult) than 
that of the strictly one-dimensional problem. Therefore, we aim at an analytical 
treatment. The electrons interact via a short-range interaction. The realistic 
case of the Coulomb interaction is supposed to be included in our treatment 
because screening removes the long-range part of that interaction. We use an 
analytical method developed recently for the study of the fc-body embedded en- 
sembles of random matrices (|ll| . Motivated by the numerical results of Ref. |^ , 
we reduce Hilbert space and keep only states where the two electrons stay rea- 
sonably close together. This leads to a non-linear sigma model. We determine 
the saddle-point manifold and thereby the influence of the two-body interac- 
tion on localization. We investigate all terms in the loop expansion up to second 
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order in the hopping matrix elements and show that these provide the expected 
renormalization of the saddle-point values. We derive a criterion for the onset 
of interaction-induced changes in the localization length. 

A brief account of this work omitting most technical details was published 
in Ref. 0. 



2 Model 

We treat the quasi one-dimensional wire in the manner introduced in Ref. . 
The wire of length L is considered as being divided into K slices labelled 
a,b,c,...,K. The surface connecting neighboring slices is transverse to the 
direction of the current through the wire. We eventually consider the limit in 
which the longitudinal length of the slices tends to zero and, simultaneously, 
K tends to infinity. In each slice the actual disorder is replaced by a random 
single-particle Hamiltonian. To simplify the notation and algebra, we consider 
the case of unitary symmetry. Hopping matrix elements connect neighboring 
slices and allow an electron to move from slice a to slices a ± 1. This model 
has been successfully used to calculate universal conductance fluctuations p2| 
and, later, the localization properties of both, the average conductance and the 
variance of the conductance |l^ . The interaction between the two electrons 
moving in the disordered wire vanishes unless both electrons occupy the same 
slice. 

In slice a we take an arbitrary basis of single-particle states labelled \aj), 
with j = 1,2, . . . ,1. We later take the limit I ^ oo. The associated creation and 
annihilation operators are denoted by al^j and aaj, respectively. The Hilbert 
space of the two-electron problem is spanned by the orthonormal states \aibj) 
defined by 

\aibj) = ala^lO) , (2) 

where |0) denotes the vacuum. This definition is unique if we require a < b and, 
for a = b, i < j . For a < b fixed, the number Nab of states is P while for a = b, 
we have Naa = ~ We always take I 3> 1 and use Naa = 

The Hamiltonian H is the sum of three terms, 

H = Ho+Hi+H2 . (3) 

The single-particle Hamiltonian 

i/o^J^/iJ^aLa,, (4) 

aij 

describes the motion of the electron within each slice under the influence of the 
disorder potential. The Hermitean matrices /i*^"-* are members of the Gaussian 
unitary ensemble of random matrices. They have dimension I, mean value zero 
and second moments given by 

h^h^;}* = ^Saa'Su'S.r ■ (5) 
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The overbar denotes the ensemble average. The parameter A has the dimension 
of an energy and is independent of the shce label a. The single-particle spectrum 
in each slice has the shape of a semicircle with radius 2A. Hopping between 
neighboring slices is described by the Hamiltonian 

aij 

It was shown in Ref. | [T^ that the hopping matrix elements u^"-* may either be 
chosen as Gaussian-distributed complex random variables with mean value zero 

and a second moment ''^ij^'^l^J* — ^Saa'^u'Sjj' or, equivalently, as elements of 
a constant diagonal matrix, 

vl^^=vS,,. (7) 

Here, v is real and independent of the slice label a. In the present context, 
we adopt the second alternative. The two-body interaction H2 is fixed (not 
random) and has the same matrix elements within each slice, 

a i<jA'<j' 

The two-body matrix elements are antisymmetric in the pairs (ij) and {i'j') 
and Hermitean. 

We wish to determine how the two-body interaction (^) affects localization 
properties. The effect (if any) will be most dramatic if both electrons stay close 
to each other as they traverse the wire. In order to focus on this point, and in 
order to simplify the calculation, we reduce the Hilbert space of two-particle 
states. We admit only states where both electrons are in the same slice, or in 
neighboring slices. Thus, of all the states introduced in Eq. (|^), we keep only the 
ones with b = a and with b = {a + 1). Without this omission, our model would 
be very general. It is, therefore, important to address the physical significance 
of this simplification. The numerical work of Ref. ||] in one dimension has 
clearly demonstrated that the change of the localization length caused by the 
two-body interaction is biggest when the two electrons which move through 
the one-dimensional wire, keep the shortest distance from each other. This is 
our motivation for the constraint artificially imposed in our model. It might 
have been desirable to loosen the constraint and to allow the electrons to keep 
a maximum distance given by the localization length in the absence of disorder. 
However, this was technically impossible. Therefore, we cannot claim that our 
results are quantitatively correct. However, in view of the work of Ref. we 
believe that we obtain qualitatively correct answers. The shortcoming of our 
approach - the reduction of Hilbert space - must be weighed against the fact 
that with this simplification, we are able to obtain a complete analytical solution 
of and a physically transparent answer to the problem. 

We aim at a comparison of the localization properties of this simplified sys- 
tem when the two-body interaction (^) is either turned off or fully present. 
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Localization properties can be read off the two-point correlation function 



C{n) = \{a{a + l)^if,\{E+ - H)-^{a + n){a + n+l)uo)\^ . (9) 

The indices /io and vq are defined in Eq. (^0|) below. To be independent of edge 
effects, we choose \ <^a<^{a + n)<^K. For large n, C{n) should decay 
exponentially with n. In proper units, the coefficient in the exponent defines 
the localization length. We calculate C(n) using supersymmetry [|l5[ 

3 Second Moment of the One— Body 

Hamiltonian Hq. Eigenvalue Expansion 

The matrix elements of the single-particle Hamiltonian in the Hilbert space 
of two-electron states (||) with b — a and b — (a + 1) are Gaussian-distributed 
random variables with mean value zero. Therefore, the distribution of these ma- 
trix elements is fully determined by their second moments. We follow Refs. |0 
and consider these second moments as matrices in the product space of two 
Hilbert-space vectors. We use this representation to derive the eigenvalue ex- 
pansion of the second moments. This will allow us to apply the supersymmetry 
technique. 

The operator i?o does not change the number of particles per slice. There- 
fore, the operator al^/^ctak' occurs in the following three types of matrix elements: 
{{a - l)iaj\Ho\{a - l)i'aj'), {aiaj\HQ\ai'aj'), and {ai{a + l)j\Ho\ai'{a + 1)/), 
with a = 1, . . . , K. Here and in the sequel, we write the adjoint of the Hilbert 
vector \aibj) as {aibj\. We simplify the notation by grouping the state labels 
into pairs, 

M = («iJi), cr = p = (^2^2), V = (i^j^) . (10) 

Thus, the Hilbert vector \aibi) is written as |a6/i). In the sequel, we refer to the 
set of states |a6/i) encompassing all values of /i as to a box. 

There are nine types of terms contributing to the second moment. Three of 
them have "diagonal form" and are given by 

—T ({a - l)an\Ho\{a - l)acr)({a ~ l)ap\Ho\(a - l)ai^) 
= ■^i^'^lla = '^{ia- l)afi\al^ack'\{a - l)aa) 

ckk' 

x{{a ~ l)ap\al,^,ack\{a - l)aiy) , (11) 

{aafi\HQ\aa(7){aap\Ho\aaiy) = ^[fj.p^ 
= ^{aafi\alf.aak'\aaa){aap\al^,aak\aai^) , (12) 

kk' 
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— {a{a + l)n\Ho\a{a + l)a){a{a + l)p\Ho\a{a + 

^■^^i^^lla = ^(a(a + l)/i|a^;,acfe'|a(a + l)cr) 

ckk' 

x{a{a+ l)p\alyack\a{a+ 1)1^) . (13) 

The sum over c in Eqs. ( |Tl| , p^ ) runs over ((a— 1), a) and (a, (a+ 1)), respectively. 
We observe that A^'i^^pa- is obtained from ^|,^^jo- by shifting a — > (a+ 1). There- 
fore, we consider only A^^i^la- and A'"^}-pa in the sequel. The "non-diagonal" 
contributions to the second moment are 

—7T {{a — l)a/i|_ffo|(a — l)a<7){aap\Ho\aah') 

x{aap\alf,,aak\aai^) , (14) 

— ((a - l)ap\Ho\(a - l)aa){a(a + l)p\H()\a{a + 
A^ 

kk' 

x{a{a + l)p\ali^,aak\a{a + , (15) 
— {aa^\Ho\aaa){a{a + l)p\HQ\a{a + 

kk' 

x{a{a+l)p\al^,aak\aia + l)i^) , (16) 

and the three Hermitean conjugate forms. 

Eqs. ( |l^ to ([l6| ) show that the matrix elements of taken in different 
boxes, are correlated. This is a trivial consequence of the fact that one of the 
two electrons occupies the same slice. The correlation causes our problem to 
differ from that of the motion of a single elctron through the wire. We return 
to this point in Sections ^ and |o| below. 

In order to be able to apply the supersymmetry technique, we now construct 
the eigenvalue expansions of the matrices A introduced in Eqs. (O) to (p^. We 
essentially follow Ref. Q. We explicitly consider Al^Ha and Ap^pa^^ and only 
state the results for the remaining matrices A. 

The matrix A\ii}},a is Hermitean, 

C /i(-i) Y ^ il(-i) (IT) 
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Therefore, can be expanded in terms of its eigenvalues and eigenvectors. 

The eigenvalue equation reads 

E4""-;p- = A(^)(-l) . (18) 

per 

Here s labels the different eigenvalues, and r distinguishes degenerate eigenvec- 
tors. The latter are normalized according to 

EC(r)(-l)(C(t'^')(-l))t = Ar(,_i), S^s'Srr' (19) 

where N(^a-i)a labels the dimension of the matrix space, see Section ||. If the 
eigenvectors form a complete set, the eigenvalue expansion of A^"^"^ is given by 

^(a-l)a 

According to Ref. we expect that the sum over s extends over s = 0, 1 only. 
This is because Hq is a one-body operator and we deal with two electrons. (It 
was shown in Ref. |11 that for a general fc-body operator and m Fermions, the 
sum over s runs from to (m — fc)). Moreover, the eigenvectors to s = are 
expected not to be degenerate and generically to have the form cjS oc S^i, = 
^iii^^hii ^hi^^ eigenvectors to s = 1 are generically expected to be (/^ — 1)- 
fold degenerate and to be given by c]mI\~\) oc ((a — l)apL\a\^acm'\{o, — \)av), 
with c = (a — l),a. The label r stands for all combinations of indices {m^m!) 
with m ^ m'. For m = m' , we have Cj^j\ — 1) oc {{a — l)aii\{al^acm — 
(1/0 X)n=i '^cmQ^cm) I (a " l)ai/) with c = (a — l),a. Using this ansatz, we do 
indeed solve the eigenvalue equation. The eigenvalues are 

A(°)(-l) = 21 ; A(i)(-1) = {21 - 1) . (21) 

Eigenvalues with s > 1 vanish. Counting shows that the eigenvectors form a 
complete set. Thus, the expansion (|o|) is established. 

We turn to A^u^lp'^ in Eq. (14). This matrix does not have the property (p^. 
Nevertheless, an expansion analogous to Eq. ( pO| ) can be found. We first deter- 
mine the right-hand and left-hand eigenvectors. It is obvious that eigenvectors 
to s = do not exist. The right-hand eigenvectors C(-^'^^(— 1, 0; r) to s = 1 are 
proportional to ((a — l)afi\a^^^_^-^^aam'\o-ai^) and belong to eigenvalue 

A(i)(-1,0) = ? - 1 . (22) 

The left-hand eigenvectors C^^^\-l,0;l) to the same eigenvalue are propor- 
tional to (aa/i|a]j„ja(a-i)m'|(a — l)ai^). All eigenvalues with s > 1 vanish. Each 
pair of left- and right-hand eigenvectors fulfills orthogonality relations analo- 
gous to Eq. ( p^ ) and can be normalized correspondingly. As a result, 
can be written in the form 

Ai-lp"J = ^(_\^ 0) - 1) E 0; i)cg^H-i. O; r) . (23) 
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The factor N is the normalization factor. For the remaining five matrices A, we 
proceed analogously. The eigenvectors and eigenvalues are listed in Appendix 
1. In summary, all matrices .4^*' with 

I = -1, 0, +1, (-1, 0), (-1, +1), (0, +1), (+1, 0), (+1, -1), (0, -1) (24) 

possess the expansion 

4^^;p. = T^EA^'^«Ec'^t-^^(*;^)^p'-'H*;r-)) . (25) 

S=0 T 

We observe that the left-and right-hand eigenvectors are Hermitean adjoints of 
each other. In the sequel, we will need the explicit values of the normalization 
factors iV(0) = f/2 and 7V(1) = f. 

4 Supersymmetry. Generating Function 

A discussion of the localization properties of the two-electron system requires 
the knowledge of the average two-point function. We calculate this function 
using Efetov's supersymmetry technique and the notation of Rcf. . The 
technique has found wide applications in the theory of disordered solids. There- 
fore, we do not give many details. We confine ourselves to indicating where the 
present treatment differs from that of Refs. 1TT| . 

The generating functional Z is an integral over commuting and anticom- 
muting variables. These are arranged in the form of the graded vectors 'i'abfi;a 
which carry the labels (ab) of the boxes, the running label fi within each box, 
and the label a which for every value of (a6/i) distinguishes two complex com- 
muting (a = 1,3) and two anticommuting (a = 2,4) integration variables. The 
functional Z contains in the exponent of the integrand a term arising from Hq. 
This term has the form 

- ^ E E Kb,L'^'{'^b^,\Ho\ab,.)L'/H,,,, . (26) 

ab A'f 

The diagonal graded matrix L is given by diag(l, 1,-1,-1). The graded indices 
have been omitted in Eq. (p6|). After averaging over the ensemble {Hq}, the 
generating functional contains in the exponent a term given by 1/2 times the 
mean value of the square of the expression (^). This term is worked out in 
Appendix 2. As a result, we find 



E E Kb,L'/'{^b^,\Ho\ab^.)Ly^^ab^ 

ab l-iv ' 

= 4 E trg[(A(-) (0))^] + ^ E trg[(A(''('^+i)) (0))^] 
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K K-l 



a—l T a—1 T 

a=2 T 

1 

a=2 T 

+ - ^ 7V(l,0)-i^trg[A(''('^+i)'")(0,l;0^^"*"+'^'"H0,l;O • (27) 



The symbols A appearing on the right-hand side denote graded matrices which 
are defined in Appendix 2. Each of the matrices A is bihnear in the integration 
variables As a consequence, the expression (|2^) is quartic in the \l/'s. To 
obtain bilinear expressions and work out the ^'-integrals, one usually performs 
a Hubbard-Stratonovich (HS) transformation. Here, this step is slightly more 
complicated than is the case usually because some of the terms in expression ( p7[ ) 
contain bilinear expressions in (rather than squares of) the matrices A. To 
overcome this problem, we recall that - with proper identification of the indices 
T - the matrices C„^ (— 1, 0; r) and Cj^l{~l, 0; I) are Hermitean adjoints of each 
other, see Eqs. (101). Therefore, the matrices 

^((a-i)a;r) ^ 1 [^(("-i)'''^) (-1, 0; r) - A^^-^)"'^) (-1, 0; /)] , 

^{{a-l)a-r) ^ 1 [^((a-l)a;r) q. ^) _^ ^((a-l)a;r) g. (38) 

2i 

have the same symmetry properties as the matrices ^^'^"^(0). On the other 
hand, using the cyclic invariance of the trace, we have 

trg[^«'^-i)"'^)(-l,0;/)A«"-i)'''^)(-l,0;r)] 

= trg[(iD«''-i)'''^))2 + (i£;(('^-i)'^'^))2] . (29) 

Therefore, the usual HS transformation can now be used for D and E. We pro- 
ceed correspondingly for A(('^-i)("+i):'")(-l, 1; r) and A^'^^'^+i):^) (0, 1; r). For 
every possible combination of slice labels {a,b) and for each of the graded ma- 
trices A we introduce a corresponding graded 4x4 cr-matrix with the following 
correspondence: 

^(ab:r) ^ ^{ab:T) 

E'^'^b-'^) ^ cr('''''^'(2) . (30) 

The volume element for integration over all these variables is denoted by d[cr]. 
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As a result, the ensemble average of the generating function Z takes the form 

Z= I d[a]exp[-^trg[^ (^(-))2 + ^ (^(a((a+l))^2 

T ^ 



2 



2 



4 ^ ^ " • 4 

i=l 



(31) 



Here N(J) is both, a graded matrix and a matrix in Hilbert space with basis 

vectors |aa/i) and |a(a + l)z/). In order to clearly display the structure of N(J), 
we list separately several types of matrix elements. There are two types of 
diagonal matrix elements, the first type given by 



(aa/x|N(J)| aav) = (e - Aa('"')^(5^ 

_^A<7(-^-)C(i-)(0) -to^. . (32) 

r 

As usual, we denote by E half the sum of the energy arguments of both Green's 
functions. The difference vanishes as both energies are equal. The matrix J 
denotes the source matrix. The second type of diagonal matrix element is 

(a(a + 1)m|N(J)| a{a + X)v) = (e - Aa("("+i»^ 5^, 

-^Xa^''^''+^^'^^Cj^j\l) . (33) 

r 

The non-diagonal matrix elements of N(J) connect either neighboring two- 
particle states, 

((a — l)a/i|N(J)|aaz/) = —v 5^^ 

-(1/2) XC^^J\-1, 0; Z)(-ic7«''-^)":^'(l) + ct((«-i)«;^)(2)) , (34) 

r 

(aa/x|N(J)|a(a + 1):^) = -v Sfj,„ 

T 

(aa//|N(J)|(a — l)azv) = —v 

-(l/2)^AC(i-)(-l,0;r)(za«»-i)«'-)(l) + a((»-i)«'-)(2)) , (36) 
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{a{a + l)fi\N{J)\aav) = -v 5^,, 



-(l/2)^AC(l,-^)(0,l;r)(z(T(''(°+i)'")(l) + (T(°(''+i)'")(2)) , (37) 



r 



or next-nearest states, 



{{a ~ l)a^\H{J)\a[a + l)u) 



= -(l/2)^AC(i-)(-l,l;0 



r 



( 



((a-l)(a+l);r)(^-^^ 



_^0-((a-l)(a+l);T)^2) 



(38) 



((a+ l)a/i|N(J)|a(a- l)t/) 



= -(l/2)^AC(i,-)(-l,l;r) 



r 



( 




_^0-((a-l)(a+l);r)(-2-) 



(39) 



The matrix N(J) also contains the source terms J. These depend upon the 
observable we are interested in. For the case of localization properties of the 
two-point function, we define the source terms in Section ^ All other matrix 
elements of N( J) vanish. 

5 Saddle— Point Equations 

For an approximate evaluation of the averaged generating function Z, we follow 
standard procedure and use the saddle-point approximation. This is suggested 
by the occurrence of large factors P and N{i) in front of the quadratic terms 
in the exponent of Eq. (|l]). We neglect J which is infinitesimal. Following the 
work of Ref. [|l2j, we also neglect the hopping terms v. For the single-electron 
problem with disorder, it was shown in Ref. that by treating w as a small 
perturbation of the saddle-point solution, the standard non-linear sigma model 
for localization is obtained. Technically speaking we assume that u <C A. The 
matrix obtained from N(J) after all these omissions is denoted by N'(0). 



In the exponent of Eq. (|3l|), we replace every cr-matrix by cr + (5cr. We collect 
the terms linear in 5a and put the coefficients multiplying the independent (5ct's 
equal to zero. As a result, we find 



The index {aa\aa) indicates that the diagonal aa box of the inverse of N'(0) has 
to be taken. The trace is taken with respect to the state indices /i, v in this box. 
Quite generally, we use the expression "box" to denote the totality of Hilbert 





(40) 
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states obtained by putting the two electrons in slices a and b, respectively, with 
b = a or b — (a + 1). Using analogous notation, we find 



(aa;r) 



P 



tr 



1 



^(a(a+l)) = ±tr 



N'(0) 

1 



{aa\aa) 



Ja(a+l);r) 



2 



{a{a+l)\a{a+l)) 

1 



N'(O)/ (a(a+l)ta(a+l)) 



A 

AC(i^)(l) 



(41) 



For the terms which are not diagonal in the box indices, it is convenient to 
introduce the matrices 



5](a6;r) ^ -crC"''^^) (1) - icr(°'';^) (2) , 



(42) 



The saddle-point equations for the index combination ((a — l)a) are given by 



2i 



iV(-l,0) 
2i 

iV(-l,0) 



tr 



1 



tr 



^'(0)/ ((a-l)a|aa) 
1 

N'(0) 



AC(i^)(_l,0;r) 
AC(i^)(-l,0;0 



(aa| (a— l)a) 



(43) 



Corresponding equations are obtained for ((a — l)(a + 1)) and a{a + 1). We test 
these results by adding Eq. ( |40| ) and the first of Eqs. (^), and the second and 
thud of Eqs. (0). With the help of Eq. (Bsl), this yields 



Aa('^'^)<5^, + J2 Aa('''''")C^i/)(0) 

T 

+ ^ Aa(°°'")C^i"Hl) 

r 

Similarly, we find 



A2 



1 



N'(0) 
1 



(44) 



I 



0) 



J2 0; r)(S«°-i)'''"))t 



p„ / {{a-l)a\aa);<jp 

A^ _ 

per 



1 

N'(0) 



(45) 



((aa| (a— l)a) ;a"p 



and completely analogous relations for the index combinations ((a — l)(a + l)) 
and (a (a + 1)). Eqs. (E3) and MB) provide a test for our calculation in the 
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following sense. On the left-hand sides of these equations, there occur the same 
linear combinations of the cr-matrices as in the matrix N'(0), see Eqs. ( p2[ ) 
through (p9[). On the right-hand side, we encounter traces over the inverse of 
N'(0) times the corresponding matrix element of the second moment of Hq. 
This is exactly the same structure as encountered by Benet et al. . We note 
that using the orthonormality relations of the coefficients C'^^'^-', we can retrieve 
the original saddle-point equations [ Eqs. (^0|) through ( ^3| ) ] from Eqs. (Q) and 
(HI). Thus, the latter set of equations encapsulates the saddle-point conditions. 

The structure of the saddle-point equations is displayed most clearly when 
we introduce the following notation. For the matrix N'(0), we write 

(a6/L<|N'(0)|a'5V) = {ES^^, - w^^ 6 at) Saa'Sbb' ~ {abp\X\a'b'iy) . (46) 

By definition, the matrix X contains all the tr-matrices which occur in N'(0). 
Needless to say that the possible index combinations (a, 6) and {a',b') are 
severely restricted by our choice of admissible Hilbert vectors. Using this nota- 
tion and the definitions of the matrices A in Section ^, we can cast the totality 
of saddle-point equations in the form 



{abfilXla'b'i') = V {abn\Ho\aba)( ) {a'b' p\Ho\a'b'i^) . 

^ V^-^"-^/(ab|a'b'):'^P 

(47) 

Eq. (p7|) is the generalization of the Pastur equation to the present problem. 



6 Solution of the Saddle— Point Equations 

We follow Ref. [0 and solve the saddle-point equations by iteration. In 
Eq. (^|) , iteration generates a continued-fraction expansion which we break off 
after n steps where n is integer but arbitrary. The n*^ element of the continued- 
fraction expansion has the form 

Y,{fxbp\Ha\abo)[-^^ {a'b'p\Ho\a'b'iy) . (48) 

pa ^ / (ab\a'b');ap 

But {E — w) is diagonal in the box indices (ab) and {a'b'). Therefore, expres- 
sion (^) equals 

'^{abp\Ho\ab(T)(—^—\ {abp\HQ\abv) Saa'Sbb' ■ (49) 

pa ^ J {ab\ab)-ap 

Using this result in the [n — Vf*' element of the continued-fraction expansion, we 
obtain a denominator which once again is diagonal in the box indices {ab) and 
{a'b'). Continuing the argument all the way up to the first term, we conclude 
that all contributions to X which are not diagonal in the box indices, vanish. 
Thus, we have for i — 1,2 

^«a-l)a;r)(^) - ; CT^'^-l^^^+l)'") (i) = ; (T^IJ^'^+'^'^Hi) = . (50) 
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The index sp stands for saddle point. It follows that the saddle-point equations 
separate into a set of {2K — 1) equations for the box-diagonal elements of X. For 
(ab) — (a(a + 1)) we have w = 0, and the n}^ element of the continued-fraction 
expansion is given by 

(a(a + l)fi\Ho\a{a + l)c7)E-^S^p{aia + l)p\Ho\a{a + 

pa 

= yA(°)(l)5,. . (51) 

The right-hand side of this equation follows from Eq. (^ ) and from the fact 
that cj:^j\—l) is traceless. The fact that E is multiplied by the unit matrix 
5^1, implies the same statement for the n*'' element of the continued-fraction 
expansion. Again, we can continue the argument all the way up to the first 
term and conclude that 

^(a(a+l);r) ^ q (52) 

Therefore, the only remaining matrix is o'("(''+i)). We use Eqs. and (p^, 
the fact that in the latter only the term with s = is relevant, and the fact that 
Cjn) = 6fi^. It follows that (t'"^"+^^) obeys the standard saddle-point equation 

^s^P^"+'^' ■ (53) 



As usual, the solution is obtained in two steps. First, we determine the diagonal 



elements 0-^°^''+^^^ which obey Eq. (|53|). For \E\ < 2A, these are given by 



rf'^+l» = |±*Ai(i?) (54) 



where we have defined 



A,{E) = \ll-{^y . (55) 



We observe that the quantity Ai{E) is proportional to the average level density 
Pspi(E). Eq. ( p5| ) shows that pspi{E) has the shape of a semicircle with radius 
2A. The invariance of the saddle-point equation Eq. ( |5^ ) under pseudounitary 
graded transformations T implies that the saddle-point manifold is given by 

^(a(a+l)) _ j,{a{a+l)) ^ i Ai{E) L ] (T^"^'^+-^^^ )~-^ 

^ 2 A 

= ^-iAi{E) T('^(''+i))L(r(''("+i»)-i . (56) 

For the index combination (aa) , the saddle-point equation contains the two- 
body interaction w. We are interested in comparing the cases w = and w ^ 
and, therefore, treat both cases. For w = 0, we conclude as before that 

= for u; = . (57) 
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For cr^""), we obtain the same saddle-point equation [Eq. (p3|)] as for a^°-''°'+'^'>\ 
Eq. ( ]40| ) and the second of Eqs. (^l]) actually differ by a factor 2. The factor 
disappears, however, when we take the traces over 6pa because the dimensions 
Naa and A^a(a+i) differ by that same factor. The solution is accordingly given 

by 



(aa) rpi^aa 




for w = . (58) 



For w 7^ 0, we encounter a novel situation because w^n is a genuine matrix and 
the arguments used above do not apply. The rfi^ element of the continued- 
fraction expansion has the form 



(aa^|i/o|aocr) [ ) {aap\HQ\aav) . (59) 



{aa\aa)\a p 



pa ^ 

The second moment of Hq in expression ( ^9| ) is proportional to the matrix 
in Eq. (^ for which we use the eigenvalue expansion Eq. (p5|). We 
first consider the terms with s = 1 for which the sum over (p, ct) takes the 
form J2pa Cp^'r''\'^)i^/[E -M){aa\aa)-ap with C^i'^^ (0) oc (aa^|aj„a„ |aat/) (where 
m 7^ n) or Cl]j\o) cx (aa/z|(aJ,jQ;„i - (l/Z) "n"")!'^'^^)- We write M^^ = 
{{E - w)-^)^,. For Cil^\0) oc (aaM|aJ„a„|aa:/), we find Epa Gii"\o)M^, = 
{^/P)Ej Mnjmj where we have returned to the original notation in Eq. (^. 
Expanding M in powers of w and taking the linear term, we note that the 
matrix element Wnjmj vanishes for m ^ n. This can easily be seen by using 
in each slice periodic boundary conditions in longitudinal direction and follows 
from momentum conservation. The same statement applies to all higher mo- 
ments of w. We conclude that Ep^{aap\al^an\aaa)M„p = 0. We turn to 

Cl^j\o) cx {aap,\{al„am - (1/0 J2n (xl,an)\aaiy) and find that J^pa CpV\o)Mp^ 
is proportional to {4:/P){J2j ^mjmj — (1/0 ^^njnj)- We again consider first 
the term linear in w. Although we do not expect the condition Wmjmj — Wnjnj 
to hold for arbitrary values of m and n, we believe that averaging over j leads 
to an approximate fulfillment of the condition Wmjmj = j ^njnj for ar- 
bitrary values of m and n. This is expected to be true a fortiori for higher 
powers of w. Thus, in very good approximation the terms with s = 1 do not 
contribute to the expression (|59|), the argument applies equally to all terms in 
the continued-fraction expansion, and we obtain 

^(aa;r) Q for W 7^ . (60) 

We are left with cr'^"''^ which obeys the saddle-point equation 

= Naa' tr[(- . (61) 

E ~~ W — ACTsp 
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This is easily verified by comparing the continued-fraction expansion of Eq. ( |6ll ) 
with that of Eq. (^9|) for a — b after we drop in the latter the terms with s = 1. 
We recall that Naa = ~ The matrix w is Hermitean. We denote the 

eigenvalues of w by Sj with j — 1, . . . , Naa- Without loss of generality, we choose 
£i < £2 < ■ • • < ^Naa- The saddle-point equation Eq. ( |6ll ) takes the form 

j E - Ej - Xa^p 

As is the case for w = 0, this equation is invariant under all graded pseudouni- 
tary transformations T^°-°'\ Therefore, we first determine the scalar solutions 
a'^"'^ and from there construct the saddle-point manifold as in Eq. ( |58| ) . Eq. ( p2[ ) 
can be written as a polynomial of order Naa in ctj""' and, therefore, possesses 
Naa real or complex solutions. We show that of these, Naa — 2 are always real. 
The remaining two solutions are either real or complex conjugate to each other. 
The positive imaginary part determines the level density pspo{E) in the box (aa). 
Therefore, the spectrum extends over that energy interval where two complex 
conjugate solutions of Eq. (p2h exist. We determine the spectrum geometrically. 



As a preparatory step, we return to Eq. (53). With z — E/X and r = 
z — a°j^°'~^^\ we write this equation in the form 

z-r = l/T. (63) 

The right-hand side is a hyperbola with asymptotes along the abscissa and the 
ordinate. The left-hand side represents a bundle of straight lines which intersect 
the abscissa at an angle of — 7r/4. For sufficiently large values of \E\, each of 
these straight lines intersects the hyperbola twice while for small values of \E\, 
there are no points of intersection. The end points of the spectrum coincide 
with the points where a straight line osculates the hyperbola, i.e., where the 
derivative of the hyperbola equals —1. This is the case at r = ±1 or, with 
Eq. (|63|), where z = ±2 and, thus, E = ±2 A. This agrees with Eq. (|5^). 

We apply the same consideration to Eq.(p2) and again use z ~ E/X and 



a^""'. Then, 



- ej/X 



Each term on the right-hand side of this equation is a hyperbola with asymp- 
totes along the abscissa and along a straight line parallel to the ordinate which 
intersects the abscissa at r — ej/X. The left-hand side again represents a bun- 
dle of straight lines which intersect the abscissa at an angle of — 7r/4. Figure 1 
represents the situation schematically. We see that for sufficiently large values 
of \E\, each straight line intersects the sum of the hyperbolas exactly Naa times, 
giving rise to Naa real solutions. For sufficiently small values of \E\, the number 
of intersections is Naa — 2. The spectrum extends continuously from a point 
El to the left of ei to a point E2 to the right of En^^- Again, the energies i?i.2 
are defined in terms of ri and T2. The latter are those points where one of the 
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straight lines z — t osculates the sum of the hyperbolas, i.e., where the gradient 
of the latter is —1. With increasing strength of the two-body interaction w, the 
difference SMaa — £i grows monotonically and so does, therefore, the difference 
E2 — El. This difference gives the width of the spectrum in the box (aa). Thus, 
the Wigner semicircle appearing in Eq. (|5^) is gradually deformed by becoming 
wider and flatter (the area remains the same). 

We supplement these general considerations by calculating position and 
width of the spectrum perturbatively for small values of w. It is easily seen 
that the first-order term only shifts the center of the semicircle from zero to 
tr{'w)/Naa, without changing either radius or shape of the semicircle. We ac- 
cordingly define Eq = ti{w)/Naa and Zq = Eq/X. Introducing the normalized 
variance of w as = {(l/iVaa)tr(w2) - [(l/iV„„)tr(w)]2}/A2 and expanding the 
right-hand side of Eq. ( p4|) in powers of w, we get 

1 

z-zq-t^-^ ^ + ... . (65) 

r T-^ 

To lowest order in [/^, the points ti and are given by ri_2 = ±(1 + (3/2)L''^). 
The end points i?i_2 of the spectrum are accordingly given by 

Ei^2 - £^0 ± (2 + U'^)\ . (66) 

This result clearly shows the widening of the spectrum. 

How strong a two-body interaction is needed to yield a sizeable effect on 
the localization length? A qualitative change of the spectrum in the box (aa) 
and, thus, in localization occurs whenever E12 and, thus, Eq and C/^/A are 
of the order of A, the width of the unperturbed two-body spectrum in box 
(aa). In order to be independent of the sign of w, we focus attention on one 
contribution to C/^, i.e., the term (1/A'aa)tr(w^)/A^ which must then be of or- 
der unity. For realistic two-body interactions, the expression tr(w^) may not 
exist. Therefore, we use completeness and write (1/Afaa)tr(w^) in the form 
(aa/z|w^|aa/i)av where the index av stands for an average over the states la- 
belled ^. We obtain (aa/z|w^|aa/i)av = . To eliminate A, we use Eq. (||) 

which yields A'^ — (^ij^)^- We identify the matrix elements ftj-J^ with the 
matrix elements (az|t4mp|ai) of the impurity potential Mmp- We do so because 
upon discretization the kinetic (potential) terms of a continuum model turn into 
hopping matrix elements (local energies, respectively). We replace the ensemble 
average by a running average over the single-particle states, use completeness, 
and obtain 

{aafi\w^\aa^)i,v > (ai|(V^mp)^|ai}av • (67) 

Whenever the criterion ( |67| ) is met, the two-body interaction will change the 
localization length qualitatively. Quantitative changes set in for smaller values 
of w^, of course. Because of the normalization of the wave functions, the crite- 
rion (p^ is essentially independent of the size of the slices. This is physically 



reasonable. Moreover, the criterion (67) can easily be checked in concrete cases. 



Since the screened Coulomb interaction is more or less fixed, the criterion (t 
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effectively establishes an upper bound on the strength of the impurity potential 
for interaction effects to be relevant. The criterion (|6^) has only qualitative 
significance because it is based upon the saddle-point approximation. When we 
identify the left-hand side with the mean-square matrix element t/^ in Eq. (Q) 
and the left-hand side with the bandwidth B"^ we see the close correspndence 
between our criterion and the result Eq. We point out, however, that on the 
quantitative level, our result differs from Eq. (^. This is discussed in Sections |^ 
and 

Returning to the general case, we write cr^""'' in the form 



cr 



aaj 



a{E)±iAo{E) (68) 



where Ao{E) > ioi Ei < E < E2 is proportional to the spectral density 
pspo{E). The general form of the saddle-point solution is then 

= a{E) - iAo[E) T("")L(r(""))-i . (69) 

The essential difference between the saddle-point solutions in Eqs. (^^ and ( ^6|) 
lies in the difference between Ao(i?) and Ai(i5), i.e, in the different spectral 
densities pspf){E) and pspi{E). Only for w = do we have Ao(i?) — Ai{E). 



7 Non— Linear Sigma Model. Localization 
Properties 

We now expand the effective Lagrangean in the exponent of Eq. around the 
saddle-point solutions obtained in Section ^. To this end, we write for i = 1, 2 



a 



{aa 



a 



(a(Q+l) 
(aa:T 



a 

fjia{a+l);T 
^((a-l)a;r)(- 
^((Q-l)(a+l);T)(-^ 
Ja{a+1)-T) 



fji^a) ^rp{aa)gp{aa)^rp{aa)^-l 

^(a(a+l)) ^ ji(a{a+l) gp{a(a+l) ^rp(a(a+l))-^-l 

y(aa)(5fj(aa;T)(^y(<Mi))-l 

rpi{a-l)a) g^{{a-l)a;T) ^-^^piaa)^-! 
'p{{a~l)a)g^{{a~l)ia+l);T)^-^^rp(a{a+l))yl ^ 
rf{aa) ^^(a(a+l);T) ^^•^^j,{a(a+l))\^-l 



(70) 



The quantities SP'^""'^ and are block-diagonal in a representation of 

the supermatrices where the first (second) block of dimension two corresponds 
to the first (the second) Green's function, respectively, see Ref. In writing 
Eqs. (^0|), we have used the freedom to multiply each Sa and each 6P from the 
left (right) by the corresponding matrix T {T~^, respectively). The volume of 
integration d[<7] in Eq. (^Tj) changes into the product of the integration measures 
for all the T's, SP's, and Sa's. Substituting Eqs. ( |70| ) into Eq. ( |3l] ) and expand- 
ing the exponent in powers of the (5P's and (Su's, we generate terms of order 0, 2 
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and i > 3. (The linear terms vanish approximately because of the saddle-point 
condition, approximately because of Eq. (|60|)). In the present Section, we focus 
attention on the terms of order zero. The remaining terms give rise to the loop 
expansion. This expansion is investigated in Sections ^ and ^ below. 

In addition to the i5ct's and SP's, N(J) also contains v and J, and we must 
clarify how to handle contributions arising from these terms. We first address 
the choice of J which is needed for investigating localization. We consider the 
average two-point function C{n) defined in Eq. (||). If the system is localized, 
this quantity decays for n ^ 1 exponentially with increasing n. The exponential 
decay is governed by the localization length. In the expression (||), the label a 
need not be equal to 1 nor do we put (a + n + 1) equal to the total number 
of slices. In fact, choosing 1 ^ a <C (a + n) <C is expected to display the 
localization properties without edge effects. Instead of the Green's functions in 
expression , we might have considered another expression where the bra state 
(ket state) is chosen as {aa\ (as \{a + n){a + n)\), respectively). Localization 
properties should not depend upon this choice, and this is indeed what we shall 
find. To generate the expression (^ from the source terms, we choose the matrix 
J to have only the following non-zero matrix elements. 

(a(a + l)n\J\{a + n){a + n + = Sf^^gSi,„gkaf3ja/3 
{{a + n){a + n + l)v\J\a{a + 1)^) = Sf^^gSi^^gkapjap ■ (71) 

Here fc is a four-by-four graded diagonal matrix with diagonal matrix elements 
+1 (—1) in the Boson-Boson block (the Fermion-Fermion block, respectively), 
and j is a four-by-four graded diagonal matrix with diagonal matrix elements ji 
{32) for the retarded (advanced) Green's function, respectively. The term (||) is 
obtained by differentiating Z once with respect to both ji and j2 at ji = = j2- 
In order to illuminate the role of the hopping terms v, it is instructive to 
display the structure of the source terms for the case of the saddle-point solution, 
i.e., when all the (5cr's and (5P's are put equal to zero. We define 

^Ka+l)) ^ y(a(a+l)) ^(^"('^+1)) {J.(<^(a+1)) yl 

^ y(a(a+l)) _ ^^(a(a+l))yl ^ 

= T^""'^ ([§ - <^{E) - i^oL) - w]-'),, (T(-))-i . (72) 

We put V — and expand tr trg InN(J) in powers of J, keeping terms of 
orders 1 and 2 which we denote by (0, 1) and (0,2), respectively. We expand 
the exponential in powers of J and keep the terms of second order, given by 
(0,1)^ and by (0,2). The term (0,1) actually vanishes and the term (0,2) is 
proportional to 

trgE {a{a + l)v\J\{a + n){a + n+l)p) 

^^((a+„)(a+„+i)) ^ ^ ^ ^ l)a\J\a{a + l)/i) . (73) 
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This is the standard term used in the calculation of localization properties us- 
ing the non-linear sigma model It is easy to see that there are non- 
vanishing contributions from higher orders in v/X. Since w ^ A, such contribu- 
tions can be omitted. 

It remains to work out the terms proportional to (u/A)^. We have already 
calculated the source terms and, therefore, put J — Q. The matrix Nsp(O) has 
the following non-zero matrix elements. 



(aa/i|N,p(0)| aaiy) = (^E - Xa^^^''^, 

{a{a + l)/i|N.p(0)| a{a + = (e - Aa^'^^^+i"^ , 

((a - l)afi\Nsp{0)\aaiy) = (aa^|Nsp(0)|(a - l)ai/) = (5^^ , 
(aa/i|Nsp(0)|a(a + 1)//) = {a{a + l)/i|Nsp(0)|aai/) = 5^^ . (74) 

We use the fact that o'sp"'' and crip^"^^'*'' are solutions of the saddle-point 
Eqs. ( |6^ ) and ([s^). Moreover, both the matrix v and (Xsp^"^^'''' are proportional 
to 5^1/. This yields 

-trg trlnN,p(0) - +(i./A)2 ^[(iV,, + iV,(,+i))/2] 

a 

x[trg(a(;'')cT(;(''+i))) + trg(a(;(''+i)a«'^+i)('^+i)))] . (75) 

The right-hand side of Eq. (|7^) contains the source terms (proportional to 
powers of J) and, most importantly, the term responsible for electron transport 
through the wire. With Eqs. (^) and ( |69| ) and the explicit values of Naa and 
Na[a+i)i this last term has the form 



+ (w/A)2(3;V4)AoAi ^ trg(r("°)L(T('"'))-iT('^('^+i»L(r(°(°+i»)-i) 

^ a 

+ ^ trg(r(°(°+^))L(T("(''+^)))~^T((°+^)(°+^)'L(r((°+^^(°+^^))"^) ) .(76) 

a ^ 

We recall that we are interested in comparing electron transport through the 
wire for the cases w = and w 0. We note that in the formulation of Eq. ([76|), 
the only difference between the two cases lies in the difference of the values of 
Ao(£;). For u; = we have Ao(£;) = Ai(£;) = ^\ - {E/2XY while this equality 
is violated for w ^ 0. 

Eq. ( [76| ) establishes a non-linear sigma model for the transport of two elec- 
trons through a quasi one-dimensional wire. This model is equivalent to the 
model for the transport of a single electron through the same wire studied in 
Ref. |]l2j. To see this, we recall that the model of Ref. also considered the 
wire as divided into K' slices numbered j = 1, . . . ,K'. Wc put K' — 2K — 1 
and map the 2K — 1 boxes labelled (aa) and a{a + 1) of the two-electron prob- 
lem onto the K' slices of the single-electron problem of Ref. [E2| by putting 
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j = 2a — 1 for the boxes (aa) and j = 2a for the boxes (a(a + 1). Then, the 
terms in Eq. ( fz^ ) take the form 

A" 

+ (i;/A)2(3;V4)AoAi^ trg(T(^U(T(^'))-iT(^'+i)L(T(^'+i))-i) . (77) 

This is exactly the form of the non-hnear sigma model derived for the transport 
of a single electron through a disordered wire in Ref. A similar correspon- 
dence can be established for the source terms. The localization properties of the 
non-linear sigma model Eq. (|7^) for single electron transport were extensively 
studied in Refs. [H @. The correspondence between the models in Eqs. ( [tt]) 
and ( [76|) makes it possible to use these results for the discussion of the localiza- 
tion properties of the model of two interacting electrons in Eq. ( fz^ ) . All we have 
to do is to transcribe the results of Refs. ^ into the present framework. 

In Refs. the continuum limit was taken by letting the length of each 

slice go to zero and the number K' of slices go to infinity. The same must be 
done in the framework of Eq. ( |76| ) . We take this step in Section ^ below. At the 
moment, it suffices to observe that the localization length is directly proportional 
to the coefficient (z;/A)^(3^^/4)AoAi appearing in Eq. ([77|). This statement 
then holds likewise for Eq. ([76[). We recall that Aq and Ai are proportional to 
the spectral densities pspo{E) and pspi{E), respectively, as determined by the 
saddle-point condition. We accordingly define the saddle-point approximation 
Tgp to the transport coefficient (sometimes also referred to as the microscopic 
dimensionless conductance) as 

T,p ^ 2n p,po{E) Pspi{E) . (78) 

This definition is patterned after the general definition of transport coefficients 
for stochastic quantum problems in Ref. p^ . Transport is mediated by the 
strength of the hopping matrix elements connecting neighboring groups of 
states. The dimensionless transport coefficient is symmetric with respect to the 
interchange of the two boxes (aa) and (a (a + 1)). Transport is possible only 



if the product pspo{E) pspi{E) differs from zero. Our result (79) is subject to 
this constraint. Let AEi {AEq{w)) denote the energy interval where pspi{E) ^ 
(pspoiE) ^ 0, respectively). We obviously have AEo{w = 0)) = Ai. Typically 
we expect that for w the interval AEo{w)) is both shifted and widened in 
comparison with AEi, the widening outweighing the shift. In this case, the 
energy interval where electron transport is possible is given by Ai?o(w — 0)) = 
AEi. 

The ratio of the localization lengths ({w ^ 0) for non- vanishing two-body 
interaction and C(w — 0) for vanishing two-body interaction, is accordingly 
given by 

Cjy^ + 0) ^ Pspo(i?) , . 

Ciw^O) p.p^{E) • ^ > 

Eq. (|7^) is valid within the same energy interval as the result ( |78| ) for the 
transport coefficient. The ratio (|7^) should typically be larger (smaller) than 
unity at the edges (in the middle) of the spectrum. 
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While physically very plausible, the form of Eq. ( [Tq ) has one obvious flaw. 
The transport coefficient is proportional to the approximate level densities as 
determined by the saddle-point approximation rather than to the exact level 
densities. And we know for sure that at least pspi{E) (which has the shape 
of a semicircle) differs from the exact level density. The latter is the level 
density of two non-interacting electrons each subject to a random Hamiltonian 
and, therefore, given by the convolution of two one-body densities of semicircle 
shape. This yields a non-zero density in the interval [— 4A, 4A] the shape of 
which is intermittent between a semicircle and a Gaussian. We conclude that 
Eq. ( I?^ ) cannot possibly be the exact answer. This is why we turn now to a 
study of the loop expansion (the expansion of the integrand of the generating 
function around the saddle-point solution) . 

8 Loop Expansion: Terms of order zero in v 

We investigate the terms in the loop expansion of Z. We recall that this ex- 
pansion is generated by expanding the exponent in Eq. (^Tj) in powers of the 
6a's and 6P's, by expanding the exponential containing terms of higher order 
than the second in a Taylor series, and by carrying out the resulting Gaussian 
integrations. 



Expanding the quadratic terms in the exponent in Eq. (31) in powers of the 
(5tT's and SP's is trivial. Terms of zeroth order vanish. Terms of first order 
cancel against those stemming from the expansion of tr trg InN(J). The terms 
of second order are 



_^trg(- (^p(-))2 + L (5p(a((a+l)))2 



T ^ 

2 
i=l 

+ iSa^-(-+^y'^)(^)f)\) (80) 



Attention thus focusses on the term — tr trglnN(J) as given by Eqs. (32) to 
(^9|), with the understanding that the substitutions (^^ have been made. 

We recall that we deal with the case v ^ X. Indeed, in the calculation of 
localization properties, it is customary to consider only the zeroth order contri- 
bution in v/X to the source term. In the calculation of the loop corrections, we 
follow this usage. It turns out that for a complete understanding of localization 
properties, terms of order are also important. These are studied in Section ^. 

Thus, the loop expansion is generated by putting u = in N( J) and expand- 
ing the logarithm both in powers of J and in powers of the (5cr's and JP's. We 
denote the terms generated in this way by (p, q) where p denotes the combined 
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power of the Ju's and 5P's, and where q denotes the power of J. We obviously 
need to consider only the cases where q = 0, 1, 2. The expansion is simplified 
because of the form of Eqs. ( |70| ) and (|7^): The factors T and T^^ cancel every- 
where except in the source terms J, and the expansion proceeds effectively in 
powers of ga^P, ffd^cr, and gdT~^JT. Here the g^s are defined in Eqs. ([72|). 

We evaluate first the lowest-order terms and, later, turn to the general 
form of the loop expansion. We consider the terms (p, q) with p < 4 and 
q < 2. The term (0,0) vanishes identically. As remarked above, the term (1,0) 
approximately cancels against the linear terms stemming from the quadratic 
expressions in the exponent of the generating function. The terms (p, 1) vanish 
for p < n. Upon expansion of the exponential and calculation of the Gaussian 
integrals, all terms which are odd in the i5cr's and SP's vanish. This leaves us 
with the following combinations of brackets (up to fourth order in total) . 

(0,2) ; (2,2) ; (4,2) ; (1,2)(3,0) ; (2, 2)(2, 0) ; (0, 2)(2, 0) ; 

(0,2)(2,0)(2,0) ; (0,2)(4,0) . (81) 

In writing the terms (2,2)(2,0) and (0, 2)(2, 0)(2, 0) we imply that we also ex- 
pand the second-order terms (2,0) originating from the expansion of the log- 
arithm into a Taylor series. This procedure differs from the standard loop 
expansion where such second-order terms are added to the ones which stem 
from the quadratic terms in the exponent of Eq. (|7^). We do so for technical 
reasons. We shall see later that this step does not affect our conclusions. The 
term (0, 2) was calculated in Section ^. We focus attention on the remaining 
terms. 

It is useful to work out the Gaussian integrals first (prior to taking the 
traces in Hilbert space). The Gaussian integrals can be performed using Wick 
contraction l|l8|. We use that 



trg{SaA)trg{6aB) = -trg(AB) , 



trg{6aASaB) = -trg(A)trg(B) . (82) 

These rules apply for arbitrary graded matrices A and B and can easily be 
checked by direct calculation. The symbol 6a stands for 5P or for any of the da's, 
and N is any of the large factors multiplying the quadratic terms in Eq. ( ^o|) . 
Needless to say, the pair of Sa's must carry equal indices (a, 6; t, i). 

Using these rules it is easy to see that all terms listed in ( |8l|) vanish iden- 
tically except for the term (4,2). This is because after applying the rules (|2|) 
until all SP's and 6a's have disappeared, we are left with a product of graded 
traces each involving powers of the propagators g^"""* or 5^'^''"^^''^ But each of 
these two propagators has the form Ffj^uSafj + G^^LaiB- Here, F and G are ma- 
trices in Hilbert space, and the indices a and /3 refer to the graded space. Thus, 
the graded traces of the g's vanish, and the same statement holds for arbitrary 
powers of these propagators. For later use, we keep track of the result, 

trg[(5^'^'^^)'=]=0 = trg[(g('^('^+i»)'=] (83) 
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which applies for all non-negative integers k. 

The term (4, 2) can be worked out explicitly. A completely analogous cal- 
culation was done in Ref. and is not repeated here. One finds that the 
term (4, 2) does not vanish for I oo but is of the same order of magnitude 
as the term (0, 2) worked out in Section ^ This finding is consistent with the 
general result of Ref. which states that if the interaction has rank k and 
the number of Fermions is m then the terms in the loop expansion do not van- 
ish if 2k < m. In the present case we have k = 1 and m = 2. We conclude 
that our result Eq. ( [79| ) for the localization length may be modified by the loop 
expansion. This is why we now examine terms of arbitrary order in the loop 
expansion. We do so in order to anwer the question: Does the loop expansion 
contribute terms which depend on the distance between the box (a (a -I- 1) and 
the box ((a + n){a + n+ 1)), i.e., between the end points of the Green's function 
in Eq. (|9|)? If the answer is yes, the behavior of the localization length displayed 
in Eq. ( |79| ) will change; otherwise, the result in Eq. (fzol ) is not affected by the 
loop expansion. 

Prior to performing the Wick contractions, the terms in the loop expansion 
have the form 

(Pl,9l)(P2,92) X ... X {Pn,qn) ■ (84) 

Here ^ ^ 9i = 2 while Pi = P > 2 must be even. We have omitted a binomial 
factor which is due to the Taylor expansions of both, the logarithm and the 
exponential. We refer to expressions like (|8^) as to "terms" , and to the indi- 
vidual factors {pi,qi) as to "brackets". Each bracket with = comprises a 
sum over all boxes. Moreover, each bracket contains a trace over Hilbert space. 
Therefore, each term in the sum over boxes defines one or several closed loops in 
the sequence of boxes. For instance, the bracket (3, 0) allows for the realization 



aa ] 



J2 trtrg[5«''-^)'^) AC(i^n-l,0;0 <5I]«"-i)-^) 

a r 
t' 

X ^ A C7(i"")(-l,l;r) [5S«'^-i)("+i)^^")]t ] . (85) 

t" 

Here and in the sequel, we use the notation introduced in Eqs. (^). In the 
expression (|85|), the loop consists of the sequence ((a — l)a) (aa) {a{a + 
1)) ((a — l)a). Another possible realization of the bracket (3, 0) is given by 

J2 tr trg [ g^r^SP^""^ 5(°°),5P("") .g(""',5p(°°) ] . (86) 

a 

Here the loop runs entirely within the box (aa). In general, we speak of a 
realization of a term if, together with the sequence of brackets (p^, also the 
sequence of 6P's and Jtr's is fixed. 

The i5P's and Sa's come in two classes. The first class comprises the forms 

gpiaa) ^ gp{a{a+i)} ^ ^c'(i^)(o) ^a^""^^) , and ^C(i^)(l) (5a(''('^+i)'^) . (87) 
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These forms connect states within a box to states within the same box. We 
refer to these forms jointly by the symbol 5D where D stands for "diagonal" . 
The second class comprises the forms 

r r 

C(i^)(-l, 1; I) (88) 

r 

which connect different boxes and move upward in the sequence of boxes, and 
the forms 

^C(i^)(-l,0;r) ]t ^ ^ C^^^' (0, 1; r) [ ]t , and 

r r 

^ 1; r) [ ]t (89) 

T 

which move downward in the sequence of boxes. We refer to the first (second) 
group jointly by the symbol 511 respectively). Within each realization of 

a term of the loop expansion, the number of JD's must be even, and to each 5Ti 
there must correspond a 5Y^ of the same type. 

For each term in the loop expansion, we first study all the realizations which 
do not carry any JS's and 5E'l''s. For these realizations, the brackets (pi, 1) 
vanish because the source term J connects the boxes (a(a + l)) and ((a + n)(a + 
71+ 1)) with n 7^ 0, and there is no way of constructing a closed loop from J and 
the 5D's. We conclude that only terms of the form (pi, 0) x . . . x (pn-i, 0)(p„, 2) 
possess the realizations here considered. The brackets {pi, 0) define closed loops 
within the same box, with a final summation over all boxes. The bracket (p„, 2) 
consists of a loop with p„ — fc steps within the box (a (a + 1)), followed by the 
step (a(a + 1)) — > ((a + n)(a + n + 1)) induced by J, followed by k steps within 
the box ((a + n)(a + n + l)), followed by the step {{a + n){a + n+l)) (a(a+ 1)) 
induced by the second factor J. Here, = 0, 1, . . . ,p„. In every bracket {pi, 0) 
with i = l,...,(n — 1) we split the sum over boxes into two parts. The first part 
comprises the sum over the boxes (a(a+l)) and ((a+n)(a+ri+l)) and is denoted 
by (pi, 0)i. The sum over all remaining boxes is denoted by (p^, 0)2- The product 
(pi, 0) X. . .X (p„_i, 0) is accordingly given by a contribution of the form (pi, 0)i x 
... X (p„_i,0)i and another contribution where at least one of the brackets 
{pi,Q) carries the index 2. The Gaussian average over the second contribution 
multiplied by (pn,2) vanishes. Indeed, since each 5D is multiplied from the 



left by a factor g, application of the contraction rules (82) leaves us, after all 
Wick contractions have been performed, with a product of graded traces of the 
form ( p3| ) and, thus, yields zero. We are left with the contribution (pi,0)i x 
. . . X {pn-ii 0)i(p„, 2). This contribution will, in general, differ from zero and be 
of the same order in Ng^(^g^^i^ as the bracket (0, 2). The term (4, 0) furnishes an 
example. The chain of terms (4, 2); (4, 2) (2, 0); (6, 2); (6, 2)(2, 0); (8, 2); . . . shows 
that terms of arbitrary order in p do contribute. Needless to say, we are not able 
to work out the contributions from all these terms analytically. Performing the 
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traces in Hilbert space is one of the stumbling blocks. We observe, however, that 
these contributions have the same form for any value of n. This is because the 
propagators g, the coefficients C*^^"^^ and the weight factors N in the quadratic 
forms (^0|) are the same for all boxes ((a + n){a + n + 1)). Therefore, the 
realizations which do not carry any factors and JE^ will renormalize the 
value of the average two-point function but will not affect the dependence of 
this function on n. In other words, the ratio ( [79| ) of localization lengths remains 
unaffected by such realizations of the loop expansion. 

We turn to the realizations which carry at least one factor 51] and (JE^ each 
and show that for n > 3, the contributions from such realizations vanish. For 
any such realization, we Wick-contract pairs of JE's and JE^'s until only a single 
such pair is left, leaving the SD^s untouched. In general, there are many different 
ways of performing these Wick contractions. Therefore, every realization gives 
rise to many different expressions with a single remaining pair ST,, ST,^ each. 
Each such expression can be written in terms of the brackets introduced above. 
Because of the contraction rules (^), the bracket structure is quite different 
from that of the original term (Js^). Our claim is proved if we show that every 
sequence of brackets involving two J's and a single pair 5E, JE^ vanishes upon 
Wick contraction. 

The pair (5E, (5E^ connects two different boxes which we denote by (a, b) and 
{a',b'). For n > 3, at least one of these two boxes must differ from both, the 
box (a(a + 1)) and the box ((a + n){a + n + 1)). Without loss of generality 
we assume this to be the box (a, 6). For n > 3 both (5E and 5E^ must occur 
within the same bracket. Otherwise, it is not possible to form a closed loop. 
This bracket contains either the string of factors 

. . .(5E^ |a&/i}(a6Ml(3d"^^^^*"^^)^5i''''^ \abiy) {abiy\ST^ . . . (90) 

with k zero or positive integer, or the Hermitean adjoint of the form (po|). Ac- 
cording to the contraction rules (|8^), Wick contraction of (5E and (5E^ yields 
the factor trg[{g'^'''' SD'-"-''^)'' g'^^'']. The remaining brackets contain the J's 
quadratically, an arbitrary number of SD^s, but no further (5E's or (JE^'s. Wick- 
contracting all the (5Z)'^"''^'s and application of the result ( |8^ ) yields zero since 
the box (ab) differs from both (a(a + l)) and {{a + n){a + n+l). This establishes 
our claim. 

The difference between the realizations which do not carry any factors ST, 
and ^E^ and those which do, lies in the fact that the former possess one variant 
(the combination (pi, 0)i x ... x {pn-i, 0)i {pn, 2) where the SD^s all carry the 
box indices (a(a +1)) or {{a + n){a + n + l)) while for the latter, the occurrence 
of at least one pair {ST, ST'') necessarily causes the occurrence of SD^s carrying 
a box index different from both (a(a + 1)) and {{a + n){a + n + l)). With the help 
of Eq. (|3|) , the Gaussian average over the latter yields zero while in the former 
case, the last graded trace involves the J's and, therefore, does not vanish. 

In conclusion, we have shown that the result ( [79| ) for the ratio of localization 
lengths remains valid to all orders in the loop expansion. 
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9 Loop Expansion: Terms of Order v"^ 



At first sight, the result of Section^ is somewhat puzzhng. Indeed, it was pointed 
out at the end of Section]^ that the result for the ratio of the localization lengths, 
Eq. (|79|), cannot be exact. This provided the motivation for the investigation of 
the loop expansion in Section || to zeroth order in w. It turned out, however, that 
Eq. ( |79| ) remains valid without modification to all orders in that loop expansion. 
This points to the need to investigate other loop corrections (which are definitely 
needed to yield a physically acceptable result). Such corrections must obviously 
be of order and will now be studied. They come in two classes. There are 
corrections of order to the source terms, and there are other corrections which 
are independent of the source terms. It is easy to see that corrections in the first 
class cannot correct the result Eq. ( |79| ) in a physically acceptable way. Indeed, 
the transport coefficient in Eq. ( [ZS] ) is proportional to Aq Ai. This factor must 
be corrected to become independent of the saddle-point solution. We note that 
this factor applies to all pairs of boxes [(aa), (a(a + 1)] and [((a — l)a), (aa)], 
for a = 1, . . . , if , and so must, therefore, the correction we are looking for. But 
by definition of the hopping matrix elements, a term of order which modifies 
the source term can affect only the pair of boxes immediately adjacent to the 
end-point boxes that occur in Eq. (^). In keeping with common practice we do 
not, therefore, consider corrections in the first class any further. 

It turns out that terms in the second class do provide the required correc- 
tions. To study such terms, we return to tr trglnN(J) as defined in Eqs. ( ^2[ ) 
to (39) with the substitutions (70). To calculate terms of order which are 



independent of J, it suffices to consider block-diagonal submatrices M("'^) and 
j^^{{a-i)a) q£ '[<i[0). Here M^"") contain the boxes [aa) and (a(a + 1)) in both 
rows and columns, 



{aaii\M^'"''^\aay) 
(aa/i|M('"'V(a+ 

aav) 

(a(a + l)/x|M(°°)|a(a+l)'^) 



{E - Xoil"'^ - AT('"')5p(°°)(T(''"')~i)(5^j. 
-^AC<,i/H0)r(°°)(5cr("''^^)(7''""')~^ , 

r 

T 

^/ji(aa)j^(a(a-M);r)^2i{a{a+l))^-l 

T 

^/jn(aa)j^j^(a(a+l);T)jt^2n(a(a+l))'j-l 

{E - Acr(°("+i» 

_^ji{a{a+l)) gp{a{a+l)) ^rp{a{a+l))'^-l^g^^ 

-A J2 (^l^J'' (l)r('^("+i»5a("('^+i)'") 



^^2n(a(a+l))-j-l 



(91) 
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while M^^'^ contains the boxes ((a — l)a) and (aa) in both rows and columns 
and is not displayed explicitly. For the terms of second order in v, we have 

tr trglnN(O) « ^ tr trg [In Af(°°) + In Af ((""i)"' ] . (92) 



To justify Eq. (92), we observe that contributions which are correct in second 
order in v can be generated from the matrix N(0) by dropping in that matrix 
all elements v but two. The set of matrices generated in this way, when in- 
serted into the expression tr trg In and expanded into powers of the w's, SP's 
and 6a's, produces the same terms of second order in v as tr trg In N(0) itself. 
In those matrices where the two remaining zj's do not appear in mirror posi- 
tions with respect to the main diagonal (as they do in Eqs. (|9l|)), the terms in 
the power series will not carry all JS's and [(^EJ^'s in pairs of the same type 
and will, therefore, vanish upon Wick contraction. Those matrices which do 
have the two remaining w's appear in mirror positions with respect to the main 
diagonal will have much larger size than shown in Eqs. (^l[), containing in ad- 
dition other boxes connected with the ones shown in Eqs. (3l])) by (5E's and 
[iJEjt's. The arguments presented towards the end of Section 8P(see the discus- 
sion around expression (pO|)) apply likewise in the present case, however, and 
show that contributions from such additional (5E's and [i5S]t's vanish. We are 
left with matrices which are block-diagonal, the diagonal block shown explicitly 
in Eqs. ( pl| ) being disconnected from the remaining one(s). When we evaluate 
the expression tr trg In, the two diagonal blocks contribute additively, and only 
]\^{aa) produces contributions of second order in v. We see that Eq. ( p^ ) holds 
except for terms which vanish upon Wick contraction and in the sense that the 
terms of second order generated by expanding either side, agree. 

We focus attention on In M^°-°-^; the terms InM^'""^)"^ are treated analo- 
gously. Expanding In Af^""^ in powers of v and the (5P's and (5cr's, keeping only 
the terms of second order in v (but terms of arbitrary order in the (5P's and 
6a's), we obtain a series of terms which individually have the form 

^2 i-j. ^j.gj^(aa)^2i(aa)-j-1^2i(a(a+l))-j^(a(a+l))^2^(a(a+l))-j-l^y(aa)-jj ^ ^gg-j 

The matrix A^°-°'^ has the form (we drop the Hilbert space indices (/x, v)) 

^(aa) ^ (5^"'')(A^P('^°) +^C(l^)(0)'5cr("'^'^))'=i<7^""^ 

r 

X AC(i^) (-1, 0; /)5E('^('^+i)'^) 

r 

^^^(,a(a+l))^^^p(Q(a+l)) _^ ^ (^(Ir) ^^-)J^(a(a+l);r)^fe2^(a(a+l)) 

r 

x^AC(i")(-l,0;r)[^E("('^+i);"']^.. , (94) 

r 

and correspondingly for A^°-^"-^^^\ The exponents fci, fc2, • ■ • are positive integer 
or zero. In both yl*^""^ and A'^°-'^°'^-^'>\ (5E's and (5E^'s with equal indices must 
always occur in pairs. 
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We expand the exponent of Z in powers of the terms (|9^) , keep only terms 
of second order in w, perform the Wick contractions, and re-exponentiate the 
result. This procedure is correct to order . It is easy to see that after Wick 



contraction, each of the terms (93) either vanishes or is proportional to 



{vjXf trg(r("°)L(r(''''))-iT(''(''+i»L(T("(''+i»)-i) , (95) 

and that the dimensionless constant of proportionality does not depend on the 
box labels (aa) and (a(a+l)). We conclude that our procedure does generate the 
expected renormalization of the transport coefficient (w/A)^ Aq Ai appearing 
in Eq. ([t^). It remains to determine the magnitude of the renormalization. 

Unfortunately, this cannot be done using the loop expansion. Even the ana- 
lytical calculation of individual terms beyond the ones of lowest order (let alone 
that of the entire series) is a prohibitive task. The difficulty lies in performing 
the traces in Hilbert space. Therefore, we resort to a different approach. We use 
a perturbative expansion which, in a similar context, was developed in Ref. | p7t 
and which is akin to diagrammatic impurity perturbation theory. This approach 
cannot give us an explicit expression either for the renormalized transport coefh- 
cient. It does allow us, however, to understand the physical origin and meaning 
of the renormalization terms. 

Starting point is the observation that in Eqs. (pi), only two boxes play a 
role. Instead of the original Hamiltonian H in Eq. (B), we, therefore, consider 
the projection PHP of H onto the part of Hilbert space spanned by the states 
\aafi) and |a(a + We consider the average two-point function 



\{aafio\iE+ - PHP)-^a{a + l)uo)\'' . (96) 

To calculate this expression, we first imagine using the technique developed in 
the present paper. The resulting averaged generating function would be written 
as an integral over the variables T^""' , r(a(«+i)) , JpC-^") , jp(a(a+i)) ^ g^iaa-r) ^ 
^^{a{a+i);T) j^{a{a+i)-T) ^g^{a{a+i)-T)^^ ^ Exccpt for a sourcc term which 
reflects the choice ( p6|) of observable and is irrelevant for what follows, the ex- 
ponent of the generating function would contain the term — tr trg InM^"") with 
]\/j(aa) (jgfjQgfj Eq (|9^). This fact then connects our renormalization prob- 
lem to the two-point function defined in (^6|). Put differently, determining our 
renormalized transport coefficient is tantamount to determining the transport 
coefficient connecting the boxes (aa) and (a(a + 1)) for the average two-point 
function (p6|). 

An alternative way of doing the calculation consists in expanding (E^ — 
PHP)~^ in Eq. ( p6| ) in powers of the stochastic variables in Eq. (^. Using 

the fact that the h'ff^^a are Gaussian distributed random variables with zero 

''J 

mean value, the terms in the resulting series can be averaged using Wick con- 
traction. The resulting contraction patterns come in two classes: (i) those where 
all contraction lines connect partners within the same box and (ii) those where 
at least one contraction line connects a partner in box (aa) with a partner in 
box (a(a -|- 1)). Eq. ( [l6| ) shows that the second class is not empty. Iteration 
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of the Pastur equation, Eq. (p7|), and comparison with our perturbation expan- 
sion shows that the saddle-point solution comprises all contraction patterns in 
class (i) with non-intersecting contraction lines. This is a well-known result. 
In the terminology of Section ^, the sum of the remaining contraction patterns 
corresponds to the sum of terms in the loop expansion. Realizations of terms 
in the loop expansion which are of order zero in 5H and [ST]'' correspond to 
the remaining contraction patterns in class (i). The rest corresponds to class 
(ii). It was shown in a different context that the remaining contraction 
patterns in class (i) are needed to change the average level density as given by 
the saddle-point condition, into the true level density. It is natural to assume 
that the same statement applies in the present context. Contraction patterns in 
class (ii) are not encountered frquently. They arise here because the interactions 
in boxes (aa) and (a(a + 1)) are correlated, and so must be the level densities 
po{E) and pi{E). 

As a result of these heuristic arguments, we are led to the conclusion that the 
loop expansion renormalizes the transport coefficient. The saddle-point version 
Tsp given in Eq. (^8|) is changed into 



The brackets indicate the average over the product of correlated densities. The 
ratio d79) of localization lengths changes accordingly into 



We now take the limit in which the longitudinal thickness d of the slices tends 
to zero while their number K goes to infinity. This yields 



We expect that the ratio on the right-hand side of Eq. ( |99| ) changes very 
smoothly with d and that a good numerical approximation can be obtained 
for finite slide thickness. 

While we have reached the end of what seems possible by analytical means, 
the arguments presented above suggest a numerical approach to the problem 
which is feasible and which can be used to determine both T and the ratio 
(^{w 7^ 0)/C(w = 0) for any given two-body interaction w. For each of the two 
boxes (aa) and (a{a + 1)), we normalize the level density so that J dEp{E) = 
1. We use p{E) = (i/(27r))(l/iV)tr[G+(i;) - G-{E)] where N stands for the 
dimension of Hilbert space in either box. Then we can write {po{E)pi{E)) in 
the form 



T = 2nv^ {po{E)pi{E)) . 



(97) 



C(^ ^ 0) _ {po{E)p^{E)) 



(98) 




(99) 



{po{E)p^{E)) 



1 



ATlNaaNa(a+l) 

y.{ir[GUE)~G-,[E)]tr[G\ 



a(a+l) 



{E) - G, 



a(a+l) 



m) -(100) 
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The ui-dependence rests in [G'+j(£') — ^"^(i?)]. The numerical calculation would 
draw the one-body Hamiltonians for the two particles in the slices a and (a + 1) 
at random, determine the Green's functions G'^^{E) and (i?) by diago- 

nalization of the resulting Hamiltonians in the two boxes, form the expression 
tr[G+ (i;) - G-,{E)] tr[G+^^,)(i;) - G-^+,)(i;)], and repeat the procedure 
many times so as to generate a meaningful ensemble average. In contrast to 
a numerical simulation of the full problem involving a large number of boxes, 
this approach is quite feasible. The simplification which we have achieved by 
analytical means consists in the restriction of the calculation to two boxes. 



10 Discussion and Summary 

We have studied the localization properties of two interacting electrons in a dis- 
order potential. Both electrons move within a quasi one-dimensional wire. This 
wire is thought of as being divided into K slices, with the surfaces separating 
neighboring slices perpendicular to the axis of the wire. The electrons interact 
if in the same slice. Hopping matrix elements allow each electron to move from 
any slice to either of the two neighboring slices. To simplify the problem, we 
have admitted only those states in Hilbert space where both electrons either 
occupy the same or two neighboring slices. We have argued that this simplifica- 
tion should be physically irrelevant: The influence of the two-body interaction 
on localization properties should not depend qualitatively on the omission of 
states in Hilbert space where the two electrons are two or more slices apart. 
The limit in which the longitudinal thickness of each slice tends to zero and the 
number K of slices tends to infinity is considered but turns out not to affect our 
results in an essential way. 

The analytical treatment of this problem becomes possible thanks to the 
eigenvalue decomposition of the second moment of the matrix elements of the 



disorder Hamiltonian in Eqs. (25). This is the novel technical feature of our 
approach. It allows us to use the supersymmetry technique. We calculate the 
average of the generating function, use the Hubbard-Stratonovich transforma- 
tion, and determine the saddle-point solutions in more or less standard fashion. 
The result (^8|) for the transport coefficient and for the ratio (^9|) for the local- 
ization lengths in the presence and in the absence of the two-body interaction 
w is intuitively convincing and demonstrates the influence of w. 

For the shape of the spectral density of two non-interacting electrons, the 
saddle-point solution yields the semicircle, a result which is manifestly not ex- 
act. This led us to consider loop corrections to the saddle-point solution. Herein 
lies the second technically novel aspect of our work. We succeeded in exhibiting 
general properties of all terms in the loop expansion up to arbitrary order. We 
could show that the loop corrections to the source terms do not alter the results 
for the transport coefficient and, thus, for the localization length obtained in the 
framework of the saddle-point approximation. We studied further loop correc- 
tions which arc independent of the source terms. For these terms we could show 
to arbitrary order in the loop expansion that each of the contributions which 
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are of second order in v has precisely the form needed to yield a renormalization 
correction to the transport coefficient. We were not able to calculate the coeffi- 
cients multiplying these terms and, thus, the magnitude of the renormalization 
effect. Instead, we used heuristic arguments to show that the renormalization 
correction does correspond to physical expectations. 

Our main result is embodied in Eq. (p9|). This is the third novel aspect 
of our work. We have shown analytically that the two-body interaction does 
affect the localization length. It does so because it alters the level density for 
those states in Hilbert space where the two electrons occupy the same slice. The 
result may be an increase or decrease of the localization length, depending both 
on properties of the two-body interaction and on the location of the energy E 
in the spectrum. A sizeable change occurs if the two-body interaction meets 
the criterion (|67|). Effectively, this puts an upper bound on the strength of the 
impurity potential which can easily be checked in practice. 

We draw attention to a characteristic difference between the one-electron 
and the two-electron problem. In the former case, the dimensionless transport 
coefficient is given in the saddle-point approximation by 27tv'^{Ai(E))'^ . The 
energy E is that of the electron. The result is exact because loop corrections 
vanish in the limit N — > oo. In the case of two electrons without interaction, the 
dimensionless transport coefficient has, in the saddle-point approximation, the 
form 27rw^(Ai(i?))^. The apparent similarity of both expressions is deceiving 
because now E is the total energy of the two electrons. Moreover, the result is 
not exact but modified by contributions from the loop expansion which changes 
the form of the level density into a convolution of two semicircles. We might, 
of course, specify the energy of each of the two electrons. But this is not a 
meaningful thing to do if we wish to compare the localization lengths without 
and with interaction. In the latter case, the total energy is the only constant of 
motion. 

The supersymmetry method does not apply to one-dimensional problems. 
Therefore, we cannot compare our result with what has been found numerically 
in one-dimcnsional systems. It is possible, however, to compare our work with 
the result Eq. (|l|) which is not restricted to one dimension. This expression 
does not contain the hopping matrix element v. Thus, it differs from our ex- 
pression (p7| ) in which we retain the structure typical for the Thouless block 
scaling argument with hopping between boxes but modify the level densities. 
Moreover, the arguments presented in the previous paragraph show that it is 
not straightforwardly possible to compare the localization length of the two- 
electron problem with that of the one-electron problem. We note, however, 
that we predict a change of localization length which depends upon all moments 
tr(?z;), tr(u;^), ... of the two-body interaction while the result Eq. involves 
only and is, thus, independent of the sign of the two-body interaction. 

To the best of our knowledge, this is the first time that a complete analytical 
treatment of the combined effects of disorder and interaction has been possible. 
We believe to have given a complete analysis of the inffuence of the two-body 
interaction on localization properties of two electrons in a quasi one-dimensional 
disorder potential. This statement is subject to the proviso that we have worked 
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in a reduced Hilbert space. 



11 Appendix 1: Eigenvectors and Eigenvalues 
of the matrices A^^^ 

The matrices and A'-~^'^'^ were treated in Section ^. The remaining 

matrices are treated correspondingly. Therefore, we simply list eigenvectors 
and left- and right-hand eigenvalues. 



z = A(°)(-l) = 2?; Cj^J {-l;r) = Cj^J ^ 5,^ ; 

A(i)(-1) ^21-1; Cl^jH-l;l) = (C^i^H-l; r))t 

oc ((a — l)afj,\aljj^acm' \ i<i — l)ai^) with c = (a — 1), a ; 
i = (0) : A(")(0) = 2(1 1); Cj^J{0;r) = Cfj{0;l) « 6,, 

AW(0) = /-2; C(i,^nO;/) = (C(l,^HO;r))t 

cx {aa^j,\al^aaTn'\aaiy) ; 
z - (+1) : AW(+1) - 2?; C(°)(+l;r) = ^^(+1;/) cx 5^, ; 

A(i)(+1) = 2/ - 1; Cll^H+^-J) = iCl'jH+l;r)y 
oc {a{a + l)^\al^j^acm' \a{a + l)^) with c = a, {a + I) ; 
i = (-l,0): A(i)(-1,0) = /-1; 

C,.l'^n-1:0;0 oc (aaM|a],„a(<j_i)„,,|(a - l)ai/) ; 

C^.i'^H-l^O;'') °c ((a - l)afi\a''^^_^^^aam'\aai') ; 
i = (0,-l); A(i)(0,-1) = /-!; 

C^.i'^HO, -1; I) oc ((a - l)a/i|a|^_^^^aam' |aat^) ; 

C/.I'^HO, -l;r) oc (aaAi|at„a(a_i)™' |(a - l)ai/) ; 
i = (0,+l): A(i)(0,+1) = ;-1; 

CMi'^HO,+l;0 Of (a(a + l)^|Q!|'^^^^^aQ„/|aaz/) ; 

C^i'^HO, +!;?') oc {aafi\ai^a(^a+i),n'\a{a + 1)^^) ; 
i = (+l,0): A(i)(+1,0) = /-!; 

C(i,^)(+l,0;r)cx (a(a + 1)m| 
z = (-!,+!) : A(i)(-1,+1) = ?; 

C(i^)(-1,+1;0 

oc {a{a + l)M|a(a+i)m"(a-i)m'l(a ^ l)ai^) ; 
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oc {{a - l)a^|Q;|^_-^^„ja(Q+i)„/|a(a + ; 
i = {+l,-l): A(i)(+1,-1) = ?; 

oc {{a - l)afi\a'i^^_^^^^a(^a+i)m'\a{a + ; 

C(i/)(+l,-l;r) 

oc {a{a + l)Ai|a[^^^)^a(a_i)„/|(a - l)aiy) . (101) 

All other eigenvalues are zero. In the first three cases, special attention must 
be paid to the cases where m — m' , see Section ||. For i — —1,0 + 1 the 
normalization factors of the Kronecker delta's are P,P/2,P, respectively. We 
take / ^ 1. Then, we have 

A^°^{i) = 2Z for i = -1,0,+1; 

A'-^^i) = / ah i except i = ±1 where A(i)(i) = 2/ . (102) 



12 Appendix 2: The second moment of expres- 
sion ( PB|) 

Because of the fact that only neighboring slices are admitted in the Hilbert 
space defined in Section S, the expression (Eq) reduces to 



2 

a—1 

. K 

E E Ka^.L^'^cia^i\H^\aav)L^'Haa. ■ (103) 

a—1 

After averaging over the ensemble {Hq\, the generating functional contains in 
the exponent a term given by 1/2 times the mean value of the square of the 
expression (102 ). We use the notation introduced in Section ^ and in Appendix 
1 and obtain for this term the following contributions. There are two "diagonal" 
terms given by 

~ 87 E E *a(a+l)M^*a(a+l)(T*a(a+l)p^*a(a+l)i"4MlV ' 
a—1 iifpcr 

E E Ka^.L^aaaKapL^aauA^;^^.p. ■ (104) 
a=l fivpa 

A diagonal term containing A^^^'^ does not appear and would be redundant. 
There are six "non-diagonal" terms given by 

X2 ^ 

-"^ZE E '^lapL'faaaKia-DpL'faia-D.A'^X'^^ , 
a=2 pi/pa 
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A2 ^ 

a=2 fj.L'pcr 
y2 K-1 

'87 X! X! *(a-l)aM^*(a-l)a<T^a(a+l)p-^*a(a+l)z.-4|«;pCT^ , 
a— 2 ^lupa 

'87 X! X! *Q(a+l)p-^*a(a+l)<T*('a-l)ap-^*(a-l)ai/-4p,Jpo^^ ' 
a— 2 ^lupa 

'87 X! X! *aaM-^*"aa**(Q+i)pi*Q(a+l)i.-4^!!;pi , 
a— 1 /^i^p sigma 

y2 K-1 

Kia+l)pL^aia+l)aKapL'i'aauAl''% . (105) 

a— 1 ^lupa 



For the sum of the terms in expressions ( |104[ ) and (105), we write somewhat 
symboHcally 

-^E E T.ll%o.Lc..^^o.%pLpp<i^pAf^^^ . (106) 

We have omitted the label indices a, (a ± 1) on the \E''s. These, however, are 
implied by the factors ^^'^ . We have likewise omitted the limits on the summa- 
tion over a. We have instead indicated the summation over the graded indices 
a and (3. We now use the decomposition in Eq. ( p5| ) for the ^^'^'s and obtain 
several contributions. The first one is due to the terms with s = and given by 

i5^trg[(A('^'^)(0))2] + ^ E trg[(A('^('^+i»(l))2] . (107) 

a— 1 a— 1 

The graded matrices A''^^'^ and A^^j^°'^'^^'^ are defined by 

l/2iT/ iTf* T 

^Wa+i) ^ zA^LV2^,(„+i),„vl;*(^^^^^^iV\ (108) 
The second one is due to the terms with s — \ and i = 0, i = 1. We find 

trg[(A(--))2] + ^ E E trg[(^('^('^-^^)-))^] . (109) 

a=l T a=l T 

The graded matrices A^^p''^\o) and A^^jf^'^'^ are defined by 

^Ka+l);r) ^ zA^LievI;„(,+i).„V&:(,+ i)^^Ly;C(l,^)(l). (110) 



,(a(a+l) _ -x rl/2,T, , iTf* T 

1/2 
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The terms with s = 1, i = (—1,0) and i = (0, —1) are given by 
1 ^ 

-^Af(-l,0)-i^trg[^«"-i)";^)(/)A«''-i)":^)(r) 

a=2 r 

^^((a-l)a,r)(^)^((a-l)a;r)(^)^ . (Ill) 

We have used that A^(— 1,0) = ^"(0,-1). The graded matrices are defined by 
A%-'^''-'^\-l,0;l) = *A^Ly„2*„„.„**^_,)„^^4/;c;.:(-l,0;/), 

^«a-i)a;x)(o^_^.;^ = iA ^ LV2^(„_i)„.„*:„^^Ly;ci-(0, -1; Z) , 

(112) 

The terms with s = 1 and i = (0, 1) and i = (1,0) yield correspondingly 

1 

- y iV(l, 0)-i y trg[^("("+i);") (0, 1; /)A("("+i):^) (0, 1; r) 

a=l r 

^^(o(a+l),r)(j^Q.;^^((a-l)a;r)^^^0.^)] , (113) 

with 

(XV 

^(a(a+i);r)(o^l.^) = iA ^ Lie*a„.„*:(„+i)^,,Ly;C7;i:;(0, 1; r) , 

4t^'^=^^(i,o;/) = iA^Ly„^*„„.„*:(„+i)^^<;c;ij:(i,o;z), 

^(a(a+i);.)(^^0.^^ = zAyLi42v|/,(,+i).„4'L^^Ly;ci:(l,0;r). 

(114) 

Finally, the terms with s = 1 and i = (—1, 1) and i = (1, —1) yield 

^ f^' iV(-l, l)-i ^ trg[A«"-i)("+i)'^)(-l, 1; Z)A«"-i)(''+i)'-)(-l, 1; r) 

a=2 r 

^^((a-l)(a+l),r)(^^ -1; Z)^((«-l)(«+l);^) (1, -1; r)] , (115) 
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with 



Ma-l){a+l)-T) 



(116) 

Eqs. (101) show that C^i^^(— 1, 0; Z) = C^i^^ (0, — 1; r) and correspondingly for 
the p airs i — (0,1), z = (1,0) and i — (—1,1),? = (1,-1). The definitions in 
Eqs. ( p2|JTT4|) and ( pi| ) imply that A((''-i)'^'^)(-l, 0; = A(('^-i)'''^)(l, 0; r) 
and similar relations for the remaining five pairs of A^s. Therefore, the sum of 
the terms in expressions (|ll^,|ll|,|ll|) simplifies to 

1 ^ 

- Y ^(-1' 0)"^ trg[A»'^-i)'''^)(-l, 0; 0v4(^"-^)'''^)(-l, 0; r) 

+ i ^ 7V(-1, l)-i ^ trg[A«'^-i)("+i)-^)(-l, 1; /)A((-i)(-+i)'-)(-l, 1; r) 

a=2 T 

1 

+ _ Y iV(l,0)-i^trg[yl('^(''+i)'^)(0,l;Z)yl('^(''+i)'^)(0,l;r-) . (117) 



The sum of the terms in expressions (107 109) and (117) is equal to 1/2 times 
the second moment of expression (p6|) . This is the result given in Eq. (E 



References 

[1] D. L. Shepelyansky, Phys. Rev. Lett. 73 (1994) 2607; see also O. N. 
Dorokhov, Zh. Exp. Teor. Fiz. 98 (1990) 646 [ Sov. Phys. JETP 71 (1990) 
360 ] 

[2] Y. Imry, Europhys. Lett. 30 (1995) 405 

[3] K. Frahm, A. Miiller-Groeling, J.-L. Pichard, and D. Weinmann, Europhys. 
Lett. 31 (1995) 169 

[4] F. von Oppen, T. Wettig, and J. Miiller, Phys. Rev. Lett. 76 (1996) 491 

[5] P. H. Song and D. Kim, Phys. Rev. B 56 (1997) 12217 

[6] F. von Oppen and T. Wettig, Europhys. Lett. 32 (1996) 741 



37 



[7] R. A. Romer and M. Schreiber, Phys. Rev. Lett. 78 (1997) 515 

[8] K. Frahm, A. Miillcr Groeling, J.-L. Richard, and D. Weinmann, Phys. 
Rev. Lett. 78 (1997) 4890 

[9] P. H. Song and F. von Oppen, Phys. Rev. 59 (1999) 46 

[10] J. Richert and H. A. Weidenmiiller, J. Phys. A: Math. Gen. (in press) 

[11] L. Benet, T. Rupp, and H. A. Wcidonmuller, Phys. Rev. Lett. 87 (2001) 
010601-1; Ann. Phys. (N.Y.) 292 (2001) 67 

[12] S. lida, H. A. Weidenmiiller, and J. Zuk, Ann. Phys. (N.Y.) 200 (1990) 
219 

[13] M. R. Zirnbaucr, Commun. Math. Phys. 141 (1991) 503 

[14] A. Mirlin, A. Miiller-Groeling, and M. R. Zirnbauer, Ann. Phys. (N.Y.) 
236 (1994) 325 

[15] K. B. Efetov, Adv. in Physics 32 (1983) 53 

[16] J. J. M. Verbaarschot, M. R. Zirnbauer, and H. A. Weidenmiiller, Phys. 
Rep. 129 (1985) 367 

[17] D. Agassi, G. Mantzouranis, and H. A. Weidenmiiller, Phys. Rep. 22 (1975) 
145 

[18] C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw-Hill, New 
York 1980 



38 




Figure 1: Graphical solution (schematic) of the saddle-point equation, Eq. (|6^. 
Some of the eigenvalues S j are shown on the abscissa. The r-dependence of the 
right-hand side of Eq. ( pi] ) is indicated by the solid lines. The two dashed 
straight lines represent the left-hand side of Eq. (Q) for two values of z. 
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